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ABSTRACT 

Recent studies suggest that the initial mass function (IMF) of the first stars 
(PopIII) was hkely to be extremely top-heavy, unlike what is observed at present. 
We propose a scenario to generate fragmentation to lower masses once the first 
massive stars have formed and derive constraints on the primordial IMF. We 
estimate the mass fraction of pair-unstable supernovae (SN^^), shown to be the 
dominant sources of the first heavy elements. These metals enrich the surround- 
ing gas up to ~ 1O~^Z0, when a transition to efficient cooling-driven fragmen- 
tation producing ^ IMq clumps occurs. We argue that the remaining fraction 
of the first stars ends up in f» IOOMq VMBHs (Very Massive Black Holes). If 
we further assume that all these VMBHs are likely to end up in the centers of 
galactic nuclei constituting the observed SMBHs, then ^ 6% of the first stars con- 
tributed to the initial metal enrichment and the IMF remained top-heavy down 
to a redshift z ~ 18. 5. At the other extreme if none of these VMBHs have as yet 
ended up in SMBHs, we expect them to be either (i) en-route towards galactic 
nuclei thereby accounting for the X-ray bright off-center sources detected locally 
by ROSAT or (ii) as the dark matter candidate composing the entire baryonic 
halos of galaxies. For case (i) we expect all but a negligible fraction of the pri- 
mordial stars to produce metals causing the transition at the maximum possible 
redshift of ~ 22.1, and for case (ii) ~ 3 x 10^^ - a very negligible fraction of the 
initial stars produce the metals and the transition redshift occurs at Zf ~ 5.4. In 
this paper we present a framework (albeit one that is not stringently constrained 
at present) that relates the first episode of star formation to the fate of their 
remnants at late times. Clearly, further progress in understanding the formation 
and fragmentation of PopIII stars within the cosmological context will provide 
tighter constraints in the future. We conclude with a discussion of several hith- 
erto unexplored implications of a high-mass dominated star formation mode in 
the early Universe. 

Subject headings: galaxies: formation - intergalactic medium - black holes - 
cosmology: theory 
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1. Introduction 

Recent studies have started to tackle the formation and collapse of the first cosmic struc- 
tures (often referred to as PopIII objects) through numerical simulations (Abel et al. 1998; 
Bromm, Coppi & Larson 1999; Abel, Bryan & Norman 2000; Bromm et al. 2001) based on 
hierarchical scenarios of structure formation. These studies have shown that gravitational 
collapse induces fragmentation of pre-galactic units with an initial baryonic mass ~ lO^M© 
into smaller clumps with a typical mass of lO^M©, which corresponds to the Jeans mass 
set by molecular hydrogen cooling. However, a considerable mass range (lO^-^M©) for the 
clumps seems plausible. 

Tracking the subsequent gravitational collapse of these metal free clumps is a very 
challenging problem as it requires the simultaneous solution of hydrodynamics and of (cooling 

lines) radiative transfer (Omukai & Nishi 1998; Nakamura & Umemura 1999; Ripamonti et 
al. 2001). Preliminary attempts and several physical arguments indicate that these clumps do 
not fragment into smaller units as the evolution progresses to higher densities. Independent 
studies (Hernandez & Ferrara 2001) comparing the observed number of metal poor stars 
with that predicted by cosmological models, also imply that the characteristic stellar mass 
sharply increases with redshift. Hence, there are grounds to believe that the first stars were 
very massive. 

The evolution of massive, metal free stars is currently subject to a rapidly growing 
number of studies (Fryer 1999; Fryer, Woosley & Heger 2001; Heger & Woosley 2001) which 
rejuvenate earher activity (Fowler & Hoyle 1964; Carr, Bond & Arnett 1984; El Eid et 
al. 1983; Pricke 1973; Fuller et al. 1986). As we discuss in detail later in the paper, stars more 
massive than about 260 Mq collapse completely to black holes, therefore not contributing 
to the metal enrichment of the surrounding gas. Similar arguments apply to stars in a lower 
mass window (30 Mq - 140 Mq), which are also expected to end their evolution as black 
holes. Hence, if supernovae from more standard progenitors (stellar masses in the range 8 
Mq - 40 Mq) are either not formed efficiently or occur in negligible numbers, it appears 
that the initial cosmic metal enrichment had to rely on the heavy element yield from the so- 
called pair-unstable supernovae (SN-j,^ ), whose explosion leaves no remnant. This conclusion 
can potentially be tested by studying peculiar elemental abundances, for example of heavy 
r-process elements (Oh et al. 2001). 

Metallicity is thought to noticeably affect the fragmentation properties of a gravitation- 
ally unstable gas. For example, Bromm et al. (2001) have shown that the evolution of a 
collapsing proto- galaxy depends strongly on the level of gas pre-enrichment. These authors 
argue that a critical metaUicity ^ 1Q~^Zq exists such that, above that value vigorous frag- 
mentation into relatively low mass (~ IOMq) clumps takes place, differently from what is 
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discussed above under metal-free conditions. 

The key question then concerns the interplay between the properties of the IMF and 
the metal enrichment of the gas. To better illustrate this, let us ideally suppose that the 
first stars are all formed - as suggested by numerical simulations - with masses above the 
SN^^ mass threshold, thus leading to the formation of black holes. Then, because metals 
are completely swallowed by the latter, the gas retains its primordial composition and star 
formation continues in the high- mass biased mode. A solution to this "star formation co- 
nundrum" must evidently exist as at present stars form with a much lower characteristic 
mass (f» IMq). 

In this paper, wc describe in detail a possible solution to this conundrum utilizing various 
constraints: metal abundance patterns in clusters, the mass density of super-massive black 
holes, off-nuclear galactic X-ray sources and, to a more speculative extent, the nature of 
a particular class (optically hidden) of gamma-ray bursts - we attempt to infer the main 
properties of the early stages of cosmic star formation. 

2. Fragmentation modes 

In this Section, we focus on the first episode of star formation and discuss the relevant 
cooling criteria and timescales in detail. We work within the paradigm of hierarchical cold 
dark matter (CDM) models for structure formation, wherein dark matter halos collapse and 
the baryons in them condense, cool, and eventually form stars. Once a halo of mass M 
collapses at redshift of Zc, the baryons are shock-heated to the virial temperature given by 

7;. = io^-VM3^/-^(i±^) K, (1) 

where Mg = M/ 10^ H'^Mq and /i is the molecular weight. If T^[^. ^ 10^ K [or equivalently 
M ^ 10^(1 + Zc)~^^'^h~^MQ], the baryonic gas cools due to the excitation of hydrogen Lya 
line. In the absence of metals, as expected for the very first episode of star formation, objects 
with Tyir ^ 10^ K can cool only through the coUisional excitation of molecular hydrogen. 
Hereafter, in this work, we refer to the former objects as Lya-cooling halos and to the latter 
as PopIII objects. 

The first stars form once the gas cools. The typical initial mass function (IMF) of this 
first generation of stars is still highly uncertain. Several authors have tackled this crucial issue 
through theoretical (Rees 1976; Rees & Ostriker 1977; Silk 1977, 1983; Haiman, Thoul & 
Loeb 1996; Uehara et al. 1996) and numerical approaches (Omukai & Nishi 1998; Nakamura 
& Umemura 1999, 2001; Abel, Bryan & Norman 2000; Bromm, Coppi & Larson 2001; 
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Omukai 2000, 2001; Ripamonti et al. 2001). Recent mult i- dimensional simulations of the 
collapse and fragmentation of primordial gas within PopIIIs, show preliminary indications 
that the IMF could be either top-heavy with typical masses of order ^ IQ'^Mq (Abel, Bryan 
& Norman 2000; Bromm, Coppi & Larson 2001) or bimodal with peaks at ^ IO^Mq and 
~ IMq (Nakamura & Umemura 2001). Note, that these numerical treatments cannot address 
the issue of the mass spectrum. Therefore, it is important to understand what physical 
processes set the scales of fragment masses and hence the stellar mass spectrum. 

It is obvious that much hinges on the physics of cooling, primarily the number of channels 
available for the gas to cool and the efficiency of the process (Rees & Ostriker 1977). In 
general, cooling is efficient when the cooling time tcooi = 3nkT/2A{n,T) is much shorter 
than the free-fall time tg = (37r/32G'p)^/^, i.e. tcooi ^fr; where n (p) is the gas number 
(mass) density and A{n,T) is the net radiative cooling rate [erg cm~^s~^]. This efficiency 
condition implies that the energy deposited by gravitational contraction cannot balance the 
radiative losses; as a consequence, temperature decreases with increasing density. Under such 
circumstances, the cloud cools and then fragments. At any given time, fragments form on a 
scale that is small enough to ensure pressure equilibrium at the corresponding temperature, 
i.e. the Jeans length scale, 

Rf^Xjoc c,ts cx (2) 

where the sound speed Cg = {RT/fiy^"^, T oc n'^~^, where 7 is the adiabatic index. Since 
Cg varies on the cooling timescale, the corresponding Rp becomes smaller as T decreases. 
Similarly, the corresponding fragment mass is the Jeans mass, 

M^ocni^^ocn^^/^+^i-"), (3) 

with rj = 2 for filaments and 1] = 3 for spherical fragments (Spitzer 1978). This hierarchical 
fragmentation process comes to an end when cooling becomes inefficient because (i) the crit- 
ical density for LTE is reached or (ii) the gas becomes optically thick to cooling radiation; 
in both cases, at that juncture tcooi ~ ts- At this stage, the temperature cannot decrease 
any further and it either remains constant (if energy deposition by gravitational contraction 
is exactly balanced by radiative losses) or increases. The necessary condition to stop frag- 
mentation and start gravitational contraction within each fragment is that the Jeans mass 
does not decrease any further, thus favoring fragmentation into sub-clumps. From eq. this 
implies the following condition: 

7 ^ (4) 
V 

which translates into 7^4/3 for a spherical fragment and 7 ,^ 1 for a filament. Thus, a 
filament is marginally stable and contracts quasi-statically when, tcooi ~ tg, and the gas 
becomes isothermal. Finally, when tcooi ^ ^fr or the fragments become optically thick to 
cooling radiation, the temperature increases as the contraction proceeds adiabatically. 



- 5 - 



2.1. Fragmentation of metal-free clouds 

In this subsection, we follow the evolution of metal-free primordial clumps during the 
fragmentation process. These results are based on the model of Omukai (2001). The gas 
within the dark matter halo is given an initial temperature of 100 K and the subsequent 
thermal and chemical evolution of the gravitationally collapsing cloud is followed numerically 
until a central proto-stellar core forms. 

The gas within the dark matter halo gets shock- heated to the virial temperature Tyij., 
which is typically ^ 100 K. However, after a short transient phase (irrelevant for the present 
analysis), the evolutionary track in the (n, T) plane shown in Fig. 1 (top curve, upper panel) 
provides a good description of the thermal evolution of the gas. The metal-free gas is 
able to cool down to temperatures of a few hundred Kelvin regardless of the initial virial 
temperature. This is the minimum temperature at which molecular line cooling becomes 
effective. 

In this analysis, external sources of heating are not included {i.e. external UV back- 
ground or CMB radiation). This is because at redshifts prior to reionization the UV back- 
ground field is relatively weak and inhomogeneous (Ciardi, Ferrara, & Abel 2000), hence we 
do not expect it to be important. Also, the CMB energy input is negligible for clouds with 
T > 100 K (or Z < IO^'^Zq) where the first stars presumably formed at redshifts z < 30. It 
might affect the temperature evolution of more metal rich systems shown in Fig. 1, but does 
not modify the conclusions drawn here. 

The gas within halos with T^[^. > lO'' K starts to cool through the hydrogen Lja line. 
It quickly reaches a temperature of ^ 8000 K and, at this stage, the fraction of molecular 
hydrogen formed is sufficient to activate H2 ro-vibrational line cooling. Thereafter, the gas 
within a Lya-cooling halo follows the same thermal evolution as the gas within PopIIIs. 
Independent of the virial temperature, the thermal evolution of the gas rapidly converges to 
the {n,T) track corresponding to Z — (the zero metaUicity track), shown in Fig. 1. The 
temperature of the gas decreases with increasing density, thus favoring fragmentation into 
sub-clumps. 

As the number density increases, it reaches the critical value ~ 10^ cm~^; and the 
corresponding Jeans mass is ^ 10^ Mq (empty dot in the upper panel of Fig. 1). The cooling 
time at this critical point becomes comparable to the free fall time as a consequence of the 
H2 levels being populated according to LTE [regime where the cooling rate A(n, T) oc n] and 
no longer according to NLTE [regime where A(n, T) oc n^]. The temperature then starts to 
rise slowly. 

At this stage, the stability of the fragments toward further increase in the density needs 
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to be investigated according to the condition eq. ^ above. The lower panel of Fig. 1 shows 
the density dependence of 7 (see eq. 0) for each metallicity track. For a metal-free gas 7 
lies in the range < 7 — 1 < 1/3, implying that further fragmentation is unlikely to occur 
unless the fragments are spherical. Although the gravitational evolution will probably favor 
a tendency towards spherical symmetry, this is not likely to occur until the central density 
has reached high 10® cm~^) values as seen from simulations {e.g. Abel et al. 2000). 

The only two deviations from the above range occur (see lower panel of Fig. 1) around 
n = 10^°cm~^ and lO^^cm"^. The former is a result of the thermal instability due to three- 
body H2 formation (Silk 1983; Haiman et al. 1996). However, this instability is quite weak 
and does not lead to fragmentation (Abel et al. 2000). The latter (at n = lO^^cm"^) is caused 
by H2 collision-induced emission. This instability is also weak and probably unimportant. 

It is important to stress that, even if the fragments are nearly spherical, fragmentation 
will be modest and is likely to result only in a low multiplicity stellar system. As the density 
increases, quasi-static contraction takes place [n = 10^°cm~^) until the fragments become 
optically thick to H2 lines, tcooi ^ ^fr, and adiabatically collapse to increasingly higher central 
densities and temperatures. At this stage, 7 > 4/3 and a central hydro-static core (stellar 
core, filled circle in Fig. 1, upper panel) is formed, with mass ~ lO^^M© (Omukai & Nishi 
1998). 

We stress again that, as long as the gas is metal-free, these sequence of events and 
conclusions hold independent of the halo virial temperature, i.e. for PopIIIs as well as for 
Lya-cooling halos. In Fig. 2 we show the evolution with temperature of the ionization fraction 
and of the fraction of molecular hydrogen for typical Lya-cooling halos of mass M = lO^M© 
at three different redshifts, z = 15, 20, 25. We find that the evolution of these fractions is 
independent of the mass (or, equivalently, T^ir) and virialization redshift. Molecule formation 
in the post-shock flow which follows the virialization of the gas within a dark matter halo 
has been recently investigated by Uehara & Inutsuka (2000). If the gas is fully ionized, 
molecules form through non-equilibrium recombination, leading to overall fractions that are 
much higher than in the expanding homogeneous universe. As a consequence of enhanced 
II2 and HD fractions, the gas rapidly cools to well below 10^ K and fragments. Molecular 
chemistry is important also for star formation in halos with T^j^ > 10^ K as assessed by Susa 
et al. (1998) and Oh & Haiman (2001). They find that initial atomic line cooling leaves a 
large residual free electron density which allows molecule formation up to a universal fraction 
of = 10~^. The newly formed molecules cool the gas further to ~ 100 K and the gas 
fragments on mass scales of a few IOOM0. 

Once the critical density for H2 {ucr ~ lO^cm"^ or lO^cm"^ for HD, if this is assumed 
to be the main coolant) has been reached, LTE conditions for the level populations disfavor 
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further cooling and fragmentation. At this stage, fragments have masses comparable to 
the Jeans mass corresponding to the point {ncr,Tcr) and virialize, following the evolution 
described above. Therefore, independent of the initial virial temperature, fragments are 
formed with typical masses ~ 10^ — lO^M©. 

Each fragment is characterized by a central core of 10~^Mq surrounded by a large en- 
velope of gravitationally unstable gas. The core grows in mass due to gas accretion from 
the envelope. The mass of the formed stars depends on the accretion rate as well as on the 
fragment mass (Larson 1999). In the absence of metals, radiation pressure cannot counter- 
balance mass accretion onto the core (Ripamonti et al. 2001, Omukai & Inutsuka 2001) and 
the mass of the resulting star is comparable to that of the original fragment, i.e. ^ IO^Mq. 
We will return to this crucial issue in Section p73[ Of course, if the parent cloud is rotating, 
angular momentum, if not dissipated, might prevent the collapse of the entire cloud. How- 
ever, for the arguments relevant to the present paper, it is sufficient that the overall efficiency 
is of a few % thus yielding a few hundreds solar mass stars out of the ~ lO^M© fragments 
predicted by Fig. 1. The question that needs to be addressed now is what mechanism can 
finally lead to the transition to a conventional mode of star formation, i.e. to a standard 
IMF, and what are the necessary conditions for this to occur. 



2.2. Fragmentation of metal-enriched clouds 

We now consider the effects of the presence of heavy elements on the fragmentation 
process. The discussion presented in this Section is based on the analysis by Omukai (2000) 
(see the original paper for details). Fig. 1, reproduced from Omukai (2000), shows the effects 
of metal enrichment on the (n, T) tracks for the same initial conditions and different values 
of the mean metallicity. 

In general, clouds with lower metallicity tend to be warmer because of their lower 
radiative cooling ability. As long as the clouds are transparent, cooling and fragmentation 
occur. Clouds with a mean metallicity Z ^ 10~^Zq follow the same evolution as that 
of the gas with primordial composition in the {n,T) plane. However, at Z ^ 10~'^Zq, B.2 
formation on grain surfaces enhances cooling at low density. Dust grains are well known to 
condense out the ejecta of SNe (see Todini & Ferrara 2001 and references therein). When 
the LTE-NLTE transition occurs for H2 , the cloud can still cool (though less efficiently) due 
to 01 line cooling. At densities > 10^ cm~^, heating due to H2 formation becomes larger 
than compressional work, i.e. 7 > 1 and the temperature starts to increase until thermal 
emission from grains due to energy transfer between gas and dust dominates the cooling. 
This occurs at a density n ^ 10^° cm~^, where the temperature drops to ~ 100 K and a 
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new fragmentation phase occurs. This shows up in the lower panel of Fig. 1 as the large 
dip in the 7 evolution. The mininiTim fragment mass is reached at the point indicated by 
the empty circle, when the density is ^ lO^^cm"^ and the corresponding Jeans mass is of 
order 10~'^Mq. Finally, as the density increases, the gas becomes opaque to dust thermal 
emission, fragmentation stops and compressional heating causes the fragments to contract 
adiabatically (7 > 4/3). Therefore, a critical metallicity of Zcr ~ 10'"^ Zq can be identified, 
which marks the transition point between metal-free and metal-rich gas evolution. 

When the metallicity is Z > IO^^Zq, at density ^ 10^ cm^'^ cooling is driven by 01, 
CI and CO line emission. When the NLTE-LTE transition for the level populations of 
CO occurs, fragmentation stops and the temperature increases due to H2 formation. The 
larger concentration of dust grains (assumed here to be proportional to the mean metallicity) 
leads to a significant thermal emission which is responsible for cooling the gas and starting 
a new phase of fragmentation. This stops when T^rom — T and thereafter the fragments 
contract quasi-statically until they become optically opaque to dust emission and adiabatic 
contraction occurs. Due to the enhanced ability to cool, fragmentation stops at lower tem- 
peratures and densities for higher metallicity clouds (see open circles in the lower panel of 
Fig. 1). However, the Jeans mass corresponding to the minimum fragmentation scale is 
always lO^^M© < Mp < lAf© for the metallicity range 10~"^ < Z/Zq < 1, several orders of 
magnitude smaller than for a cloud with no metals. 

At the onset of adiabatic contraction, when 7 becomes larger than 4/3, an initial hydro- 
static core forms (transient core) with a mass ~ 10~'^Mq regardless of the metallicity of 
the gas. This transient core is fully molecular and is absent when the gas is metal-free. 
The temperature of the transient core increases as its mass increases due to accretion of 
surrounding gas. Eventually, the temperature reaches about 2000 K, where H2 dissociation 
begins. This softens the equation of state until the dissociation is almost complete. Then, 7 
falls below 4/3 in the density range 10^^ — 10^°cm~^. Note that the thermal evolution after 
H2 dissociation (i.e., n > lO-'^^cm^^) is the same independent of the initial composition of 
the gas. After 7 once again exceeds 4/3, a hydro-static core (so-called stellar core) forms. 
Its physical characteristics are independent of metallicity (n ~ 10^^cm~^,M «i 1O~^M0; 
Omukai 2000). 



2.3. Formation of proto-stars 

Proto-stars, whose mass is initially very low (about 1O~^M0), grow in mass by accretion 
of the envelope material. The final mass of stars is, then, determined not only by the mass 
of the fragments, but also by when accretion stops. The accretion rate onto the proto-star 
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is related to the sound speed (and therefore, temperature) of the proto-stellar cloud by the 
relation, M ~ cf/G, (Stahler, Shu, & Taam 1980). 

As seen in Fig. 1, the temperature of proto-stellar clouds decreases with metallicity. 
Therefore, the mass accretion rate is higher for proto-stars formed from lower metallicity 
gas. If dust is present in the accretion flow, accretion onto the massive proto-star becomes 
increasingly difficTilt owing to the radiation pressure onto the dust. In present-day interstellar 
gas, accretion onto stars more massive than SOM© is inhibited by this mechanism (Wolfire 
& Cassinclli 1987). In gas with lower metallicity, the mass bound is expected to be higher 
because of the higher accretion rate and lower radiation force. In particular, for metal-free 
gas, this mechanism does not work. Therefore, without dust, accretion is likely to continue 
until the ambient gas supply is exhausted (Ripamonti et al. 2001, Omukai & Inutsuka 2001). 

In conclusion, the presence of metals not only enables fragmentation down to smaller 
mass scales lO'^M©, but also breaks the one-to-one correspondence between the mass of the 
formed star and that of the parent fragment by halting the accretion through radiation force 
onto the dust. 



3. The star formation conundrum 

According to the scenario proposed above, the first stars that form within PopIII and 
Lya-cooling halos out of gas of primordial composition tend to be very massive, with masses 
10^~^Mq. It is only when metals change the composition of the gas that further frag- 
mentation occurs producing stars with significantly lower masses. It is at this stage that we 
expect a transition to occur from a top-heavy IMF towards a more conventional IMF with 
a wide range of masses, such as the one observed locally (Scalo 1998, Kroupa 2001). 

A growing body of observational evidence points to an early top-heavy IMF (see Her- 
nandez & Ferrara 2000 and references therein) . Compelling arguments for an early top-heavy 
IMF can also be made from observations on various scales: (i) the early enrichment of our 
Galaxy required to solve the so-called G-dwarf problem, (ii) the abundance patterns of met- 
als in the Intracluster Medium (ICM), (in) the energetics of the ICM, (iv) the non-detection 
of metal- free stars, and (v) the over-production of low mass stars at the present epoch and 
metals at z ~ 2 — 5 for the submillimeter derived star formation histories using SCUBA 
detections for a standard IMF can be resolved with an early top-heavy IMF. 

For instance, the ICM metal abundances measured from Chandra and XMM spectral 
data are higher than expected from the enrichment by standard IMF SN yields in cluster 
galaxy members, which can be explained by a top-heavy early IMF. Furthermore, the ob- 
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served abundance anomalies {e.g. oxygen) in the ICM can be explained by an early generation 
of PopIII SNc (Locwcnstcin 2001). There is observational evidence from the abundance ratio 
patterns of [Si/Fe], [Mg/Fe], [Ca/Fe] and [Ti/Fe] in the extremely metal-poor double-lined 
spectroscopic binary CS 22876-032 in the halo of our galaxy (Norris, Beers & Ryan 2000) 
for enrichment by a massive, zero-metallicity supernova on comparison with the theoretical 
models of Woosley & Weaver (1995). These issues, highly suggestive of top-heavy early star 
formation, have recently motivated a series of numerical investigations of the nucleosynthesis 
and final fate of metal-free massive stars (Heger & Woosley 2001; Fryer, Woosley & Heger 
2001; Umeda & Nomoto 2002). In their recent paper, Heger & Woosley delineate three mass 
ranges characterized by distinct evolutionary paths: 

1. >, 26OM0: the nuclear energy release from the collapse of stars in this mass range 
is insufficient to reverse the implosion. The final result is a very massive black hole 
(VMBH) locking up all heavy elements produced. 

2. 14OM0 <, <, 260Mq: the mass regime of the pair- unstable supernovae (SN^^ ). Pre- 
collapse winds and pulsations result in little mass loss, the star implodes to a maximum 
temperature that depends on its mass and then explodes, leaving no remnant. The 
explosion expels metals into the surrounding ambient ISM. 

3. 3OM0 ^ ^ I4OM0: black hole formation is the most hkely outcome, because either 
a successful outgoing shock fails to occur or the shock is so weak that the fall-back 
converts the neutron star remnant into a black hole (Fryer 1999). 

Stars that form in the mass ranges (1) and (3) above fail to eject most of their heavy 
elements. If the first stars have masses in excess of 26OM0 (in agreement with numerical 
findings), they invariably end their lives as VMBHs (in the following we will refer to VMBHs 
as black holes of hundred solar mass or so) and do not release any of their synthesized heavy 
elements. However - as we have shown - as long as the gas remains metal-free, the subsequent 
generations of stars will continue to be top-heavy. This star formation conundrum can be 
solved only if a fraction of the first generation of massive stars have masses ^ 26OM0. Under 
such circumstances, these will explode as SN-j,^ and enrich the gas with heavy elements up 
to a mean metallicity oi Z ^ 10"'^ Zq, and as per arguments outlined in the previous section 
(see Fig. 1), thereafter causing a shift over to an IMF that is similar to the local one. In what 
follows, we further explore the implications of the above scenario and derive the conditions 
for the solution of the conundrum. 
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3.1. Abundance of VMBHs and metals 

As a consequence of the picture proposed above, VMBHs are an inevitable outcome. 
We now compute the expected mass density of metals and the mass density of remnant 
VMBHs produced in such a first episode of star formation. For a collapsed dark matter halo 
of total mass M, the associated baryonic mass is assumed to be M{ilB/^M) (where VLb/^m 
is simply the baryon fraction). Following the results of numerical simulations at different 
resolutions (Bromm et al. 2001, Abel et al. 2000) we assume that Ri 1/2 of the total available 
gas is utilized in star formation, the rest remaining in diffuse form. 

The relative mass fraction of SN^^ and VMBH progenitors are parametrized as follows: 

^77 = /77 V = /77^*, (5) 

M. = (l-/,,)-(^) = (l-A,)M., 

where M^^ (M,) is the total mass which ends up in SN^-y (VMBHs) and M* is the mass 
processed into stars. Thus, only a fraction /^^ of the formed stars can contribute to gas 
metal enrichment. The metal yields for the dominant elements have been computed using 
the results of Heger & Woosley (2001), and are plotted in Fig. 3 as a function of the mass 
of the progenitor star. The bulk of the yield is contributed by O^^ and the mass of metals 
ejected can be written as 

Mz ^ M™/2 - IOMq ^Z". , (6) 
°2OOM0' ^ ^ 

where we have taken f» 2OOM0 as a fiducial mass for SN^^ progenitors. Next, the further 
assumption is made that metals are ejected from the parent galaxy into the IGM and their 
cosmic volume filling factor is close to unity, therefore uniformly polluting the IGM. This 
comes from the results of Madau, Ferrara & Rees (MFR, 2001) and Mori, Ferrara & Madau 
(2001), where it is shown that for reasonable values of the star formation efficiency, metal 
bubbles produced by (standard SNII in) proto-galaxies which result from 2-a peaks of the 
density power spectrum at redshift 10 do overlap. Indeed, the kinetic energy released by 
SN^^ is much higher than for ordinary SNII (see Fig. 2), hence making our assumption 
even more solid. The temperature of the hot gas is expected to be somewhat higher than 
estimated by MFR: however, due to the higher redshifts and consequent stronger Inverse 
Compton cooling, ejected metals have enough time to cool by 2; ~ 3. 

As long as the average metallicity is below the critical value, i.e. Z^r = IO^^Zq, we 
argue that the IMF remains top-heavy and the redshift-dependent critical density of metals 
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contributed by all SN^^ at redsliifts > z can be computed, 

1 /-oo r 

nz(^z) = — dz' dMn{M,z')Mz (7) 

Per JMrnin{z') 

where n{M, z) is the number density of halos per unit mass predicted by the Press-Schechter 
formalism and the integration is performed from Mmin(-2) which is the minimum mass that 
can cool within a Hubble time at the specified formation redshift z, i.e. tcooi{M,z) ^ tui^z) 
(see Ciardi et al. 2000). We adopt a cosmological model with the following parameters to 
compute the abundance of halos : VLm = 0.3, = 0.7, h = 0.65, Qb = 0.047 (latest 
predictions from BBN, see Buries, NoUett & Turner 2001) and use the COBE-normalized 
power spectrum for fluctuations as described in Efstathiou, Bond & White (1992). 

From the above expression, we can estimate the transition redshift at which the mean 
metallicity is 10"^^©, 

< Ziz,) >= ^ ^ 10-Ze (8) 

for various values of the fraction /^^. The results are plotted in Fig. 4 alongwith the corre- 
sponding critical density flvMBH contributed by the VMBHs formed, which is given by, 

^VMBHizf) = ii-^^ r dz' f dMn{M, z') M,. (9) 

Pcrit "'^/ JMminiz') 

Several assumptions as to the fate of these VMBHs at late times will be examined in an 
attempt to obtain constraints on the value of /^^ from current observations in the next 
Section. 



4. Observational Constraints 

The objective now is to set some limits on /^^, the fraction of stars formed that result 
in SN-y-y providing the first metals, that are eventually responsible for the shift to a normal 
IMF. In order to compare the predicted critical density of VMBH remnants to present 
observational data we need to make further assumptions about the fate of these VMBHs at 
late times. 

The mass density contributed by remnant supermassive black holes ^%mbh can be 
estimated observationallyQ from the demography of local galaxies (Magorrian et al. 1998; 



The superscript denotes the present day value 



-13- 



Merritt & Ferrarese 2001). There are two extreme possibilities (bounding cases for the values 
of f^ry) concerning the relation between the inferred mass density of SMBHs and VMBHs: 

[A] ^vMBH{zf) — ^sMBH- VMBHs havc, during the course of galaxy mergers, in 
fact, been used to build up the SMBHs detected today; 

[B] ^vMBH{zf) = ^vMBH 7^ ^SMBH there is no relation between the detected SMBHs 
and the early VMBHs, implying that SMBHs have grown primarily by gas accretion. 

Thus, in the first scenario [A], we assume that all VMBHs are able to sink to the center 
of host systems to merge within a Hubble time, giving birth to the SMBHs in the nuclei of 
galaxies that we see today. Note, that this is in fact unlikely (as explained further below), 
nevertheless, this argument provides an extreme bound. 

It is important to point out here, that in the current proposed theoretical models for 
the mass accretion history of black holes, the mass density seen today locally in SMBHs can 
be built either predominantly by mergers or by accretion. In phenomenological models that 
attempt to tie in the current high redshift observations of the space density and luminosity 
function of quasars with that of the local space density of black holes wherein it is assumed 
that quasars are (i) powered by BHs and (ii) are optically bright for a period of f« 10^~^ 
yr (Haiman & Loeb 1998; Haehnelt, Natarajan & Rees 1998; Kauffmann & Haehnelt 2000), 
either picture (mergers/accretion) can adequately explain the mass assembly of the remnant 
SMBHs. 

This permits both scenarios [A] (wherein SMBHs grow primarily via mergers of VMBHs 
during galaxy assembly) and [B] (wherein SMBHs grow primarily by accretion and ^%mbh 
is unrelated to VLymbh) to be plausible extreme cases. Consideration of cases [A] and [B] 
will provide bounds on the value for /^^. For a given /^^, the corresponding value for the 
transition redshift can be found on the Zf curve shown in Fig 4. 

We note here that VMBHs in the mass range expected from the first episode of metal- 
free star formation (a few hundred solar masses) are likely to be ubiquitous. Using ISO 
observations, Gilmore & Unavane (1998) have convincingly argued that low mass stars (in 
fact hydrogen burning stars of any mass) do not provide a substantial portion of the dark 
matter inferred in the halos of galaxies. In combination with the mass limits obtained on 
potential dark matter candidates from the MACHO (Alcock et al. 2001) and EROS (Lasserre 
et al. 2000) experiments (micro-lensing studies) in our Galaxy, low mass objects, with masses 
under 10 are ruled out - thereby, making VMBHs IOOM0 BHs) in halos (case [B]) 
plausible dark matter candidates. 

For case [A], we can simply equate ^vmbh to the measured mass density in SMBHs 
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found from the demography of nuclei of nearby galaxies (Magorrian et al. 1998; Gebhardt 
et al. 2001) indicated by the horizontal line in Fig. 4 (top panel). The most recent value 
reported by Merritt & Ferrarese (2001) is ^smbh — ^O'^^hfls, which for h = 0.65 and 
Qb = 0.047 yields ^%mbh — 3-05 x 10~^. From the top panel of Fig. 4, we see that this 
gives /^^ ~ 0.06. The corresponding value for the transition redshift can be found on the 
Zf curve on the same plot, yielding a value Zf ~ 18.5. This constraint implies that all but 
a small fraction of the mass involved in the first episodes of star formation - approximately 
6% - went in objects outside the mass range 14OM0 < < 26OM0, yielding remnant black 
holes with mass 10 — 500 Mq. The scenario corresponding to [A] therefore implies that Zf, 
the transition redshift from top-heavy to normal IMF occurred at 18.5. 

Let us now look at the other extreme case [B], wherein the assembled SMBHs in galactic 
centers have formed primarily via accretion and are unrelated to VMBHs. This leads us to 
two further possibilities: (i) VMBHs would still be in the process of spiralling into the centers 
of galaxies due to dynamical friction but are unlikely to have reached the centers within a 
Hubble time due to the dynamical friction time-scale being long (Madau & Rees 2001), or 
(ii) VMBHs contribute the entire baryonic dark matter in galactic halos. 

Pursuing scenario (i), therefore, at best some fraction of these VMBHs (for instance, 
those in gas-rich regions of gas-rich systems) might appear as off-center accreting sources 
that show up in the hard X-ray wave-band. Such sources have indeed been detected both by 
ROSAT (Roberts & Warwick 2000) and more recently by Chandra in nearby galaxies (M84: 
Jones 2001; NGC 720: Buote 2001). 

In a survey of archival ROSAT HRI data to study the X-ray properties of the nuclei 
of 486 optically selected bright nearby galaxies (Ho, Fillipcnko & Sargent 1995), Roberts & 
Warwick (2000) found a large number of off-center X-ray sources. The X-ray sources detected 
within the optical extent of these galaxies were classified cither as nuclear or non-nuclear 
(and therefore off-center) depending on whether the source was positioned within 25 arcsec 
of the optical nucleus. They detect a nuclear source in over 70% of the galaxy sample and a 
total of a 142 off-center sources. Roberts & Warwick (2000) find that the non-nuclear sources 
follow a steep, near power-law X-ray luminosity distribution in the 10^^ — 10^° erg s~^ which 
leads to an Lx/Lb ratio of 

^ = 1.1 X 10^^ erg s-i(10^° Lq)-^ (10) 
Lb 

The median luminosity of the non-nuclear sources is found to be a factor of ~ 10 lower 
than that of the nuclear sources. However, they estimate the incidence rate of off-center 
sources with Lx > 10^^'^ erg s~^ (which corresponds to the Eddington luminosity for a 1.4 
solar mass neutron star) to be 0.7per lO^'^ galaxy. The existence of accreting VMBHs 
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might also help explain the following observation. The far-IR excess detected with DIRBE 
at 60 and 100 fim has been tentatively interpreted as an extra-galactic background with 
integrated intensity of 44 ± 9 nW m~^sr~^ in the range 60 — 100/xm (Finkbeiner, Davis & 
Schlegel 2000). The energy required to produce such a FIR background could derive from a 
highly obscured population of accreting VMBHs at moderate redshifts (for an alternative to 
this see explanation of scenario (ii) below). 

In order to estimate the mass density contributed by these VMBHs, we can compare 
the mass in VMBHs to that of the SMBH in a typical galaxy of luminosity say 10^° Lq for 
case [B]. Most of the VMBHs are likely to be en-route to galactic centers. Note however, 
that not all the inspiralling VMBHs will be accreting and be X-ray bright: many of these 
could, in fact, be low radiative efficiency ADAFs (Advection Dominated Accretion Flows), 
and therefore be too faint to be detected by ROSAT or Chandra or not accreting at all if they 
do not happen to be in gas-rich regions of the galaxy. Thus, using the following argument 
we can obtain a lower limit on the abundance of VMBHs. 

The mass of the central SMBH hosted by a galaxy with a luminosity Lgaiaxy = 10^° Lq 
(assumed to be the fiducial galaxy for these purposes) in the B-band is given by (Merritt & 
Ferrarese 2001): 

MsMBH = 10 ^Mbulge = 10 ^-C/galaxy 

where the mass-to-light ratio in the B-band is taken to be 1 Mq/Lq and the galaxy luminosity 
is essentially dominated by the bulge luminosity. Hence, the mass of the central SMBH is 
MsMBH ~ IO^Mq. Now, we use the luminosity of the ROSAT off-center sources (Roberts 
& Warwick 2000) and use the fact that they find ~ 0.7 off-center sources with luminosity 
> 10^^ '^ erg s~^ per lO^'' Lq galaxy, to estimate the mass of VMBHs in the fiducial galaxy. If 
we assume that these VMBHs have a typical mass of ~ SOOMq, their Eddington luminosity 
is, 

-^^Edd = 4.3 X 10 ( ) erg s . 

^ 3OOM0 ' ^ 

Thus, VMBHs appear to be radiating at sub-Eddington luminosities, the average rate (given 



by eq. |10]) being 3% of the Eddington value. This is consistent with the observed lumi- 
nosities of the nuclear sources in the sample, which Roberts & Warwick (2000) find to be 
radiating at severely sub-Eddington rates. 

Therefore, we expect that the typical accreting VMBH mass in a 10^° Lq galaxy to be 
about 210 M0. Little is known about the spatial distribution of such objects. How can we 
take into account the fact that only a fraction of the sources are likely to be accreting? Since 
the detected ROSAT off-center sources are within the optical radius of the galaxies, the 
number of non-accreting VMBHs in the halos of these galaxies can be large. One can obtain 
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a simple estimate of this number by pursuing the following argument. Assume that VMBHs, 

being collisionlcss particles, closely trace the dark matter (which we take to obey a Navarro, 
Frenk & White [NFW, 1997] profile) and that the ratio of virial to optical radius is ~ 15 for 
a 10^^ Lq disk galaxy (Persic, Salucci & Stel, 1996). These hypotheses require that we scale 
up the total mass in VMBHs by a factor of ~ 10, which gives Mvmbh = 2.1 x 10^ Mq. Now, 
the ratio Mvmbh/Msmbh is 2.1 x 10""^, implying that ^vmbh = 2.1 x 10~'^ilg]^^jf = 
6.4 x 10~^° which in turn gives f^^ ~ 1 (see dashed line in bottom panel of Fig. 4) and 
Zf ~ 22.1. 

Now we explore scenario (ii) of Case [B], the instance when VMBHs constitute the entire 
baryonic dark matter content of galaxy halos (but do not contribute to the disk dark matter) . 
It is important to point out here that cosmological nucleosynthesis arguments require both 
baryonic and non-baryonic dark matter (Pagel 1990). Essentially, this is due to the fact that 
the mass density contributed by baryons Jls, is well in excess of Qy, the contribution to the 
mass density by visible baryons. 

Using the luminosity dependent relation of visible to dark matter for spirals (Persic, 
Salucci & Stel, 1996) for a fiducial galaxy of 10^° Lq we find, 

Ri 0.05 



DM 

where M^ig is the visible mass and Mjjm = + ^d'm total dark matter mass, 

given as a sum of a baryonic and a non-baryonic component. The total baryonic mass is a 
fraction Qb/^m of the total mass of the system, e.g. 

M^is + M^j^j^ = ^{M^is + Mdm)- 

From the two equations above and assuming a NFW density profile, we estimate that the 
total baryonic dark matter content as, 

Mdm - (^M - ^)M,is = 2.33M^,,; 

moreover, 90% of this mass resides outside the optical radius and can be contributed by 
VMBHs: 

Mvmbh 2.1 x M,js / x 

^ = ^ = 2.1xl0^ 11 

Msmbh 10-3 xM,,, ' ^ ^ 

implying 

^VMBH = 2.1 X 10^ X nsMBH = 6.4 X 10~^ (12) 

Incorporating this constraint into Fig. 4 (bottom panel), we find f^j 3.15 x 10"^ and 
Zf ~ 5.4. Note that this scenario does not violate the constraint on the over-production of 
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background light (see Carr 1998; Bond, Carr & Hogan 1991; Wright et al. 1994). In fact, 
some fraction of the observed ncar-IR DIRBE excess could be produced by VMBHs. In a 
recent estimate of the cosmic background at 1.25 fim and 2.2 fim (corresponding to the J 
and K bands respectively) using 2MASS and the DIRBE results Cambresy et al. (2001), 
also find an excess (significantly higher than the integrated galaxy counts in the J and K 
bands) suggesting the contribution of other sources. Pop III stars and their VMBH remnants 
(accreting at very high redshifts) postulated here are likely candidates. 

According to scenario [B], the average metal abundance at rcdshift ~ 5.4 should be 
Zcr = IO^'^Zq, marking the transition from a top-heavy IMF to a standard power-law IMF. 
This does not necessarily violate the observed metal abundances in Damped Lya systems, 
Z Ri 10~^, as it is yet not clear to what extent the IGM metaUicity is spatially uniform at 
these intermediate redshifts. 

Therefore, scenario [B] is poorly constrained by the data, as the two options (i) and (ii) 
give large ranges for /^^ and Zf-. 3.15 x 10^^ < /^^ < 1 and 5.4 < Zf < 22.1. Consistently, 
scenario [A] gives a limit on /^^ and Zf of /^^ = 0.06 and Zf — 18.5. 

The actual data do not allow us at present to strongly constrain the above two quanti- 
ties, but they provide interesting bounds on the proposed scenario. These ranges also have 
imphcations for the expected detection rate of SNe beyond Zf with future instruments hke 
NGST (Marri, Ferrara & Pozzetti 2000). 

5. Discussion 

We have proposed a scenario to solve a puzzling star formation conundrum: the first 
stars are now thought to be very massive and hence to lock their nucleosynthesis products 
into a remnant (very massive) black hole. This high-mass biased star formation mode con- 
tinues as long as the gas remains metal-free. During this phase, metal enrichment can occur 
only if a fraction f^^ of the stars have mass in the window leading to pair-unstable SNe 
(14OM0 < < 26OM0) which disperse their heavy elements into the surrounding gas. 
Such metals enrich the gas up to Z^r ~ 10~^ Zq, when a transition to efficient cooling-driven 
fragmentation producing ^ IMq clumps occurs at rcdshift Zf. We argue that the remain- 
ing fraction of the first stars end up in ^ IOOMq VMBHs. By analyzing the evolutionary 
fate of such objects, we argue that [A] they could end up in the SMBHs in the centers of 
galactic nuclei, [B] (i) could be en-route to the center and hence identified with the X-ray 
bright, off-center ROSAT sources, or (ii) constitute the entire baryonic dark matter content 
of galaxy halos. These possibilities are used to obtain constraints on the two quantities: 
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/^^ ~ 0.06 and Zf ~ 18.5 for [A], and /^^ f« [10"^ - 1] and Zf f« [5.4 - 22.1] for case [B]. 
The value Z^r ~ 10~^Zq found here is admittedly somewhat uncertain. For this reason we 
have investigated how the above results might be affected by a different choice, e.g. assum- 
ing Zcr ~ 10~^Zq. Indeed, lower values of Zcr imply that less efficient metal enrichment is 
required in order to change to a more conventional star formation mode. Thus, comparable 
values for Zf are found but the corresponding values for f^^ tend to be systematically smaller, 
being /^^ 0.6% for scenario [A] and in the range 3.15 x 10~^ < Zf < 0.98 for scenario [B]. 

Several uncertainties remain in the comparison of the inferred density of VMBHs to 
local observations. For example, it is not obvious if SMBHs could at all form out of VMBHs. 
Dynamical friction can effectively drag VMBHs towards the center of the host system, at 
least within a distance of ~ 100 pc (Madau & Rees 2001). Unless most of the energy is 
radiated away in gravitational waves, it could be difficult for a VMBH cluster to coalesce 
into a single unit. Furthermore, accretion onto isolated VMBHs could be too inefficient to 
explain most of the off-center sources observed by ROSAT and CHANDRA. Higher accretion 
rates might be activated by the tidal capture/disruption of ordinary stars. Question remains 
if the frequency of such an event is in fact sufficiently high to explain the data. 

In spite of the many difficulties and uncertainties discussed above, our study represents 
a first attempt to link the first episode of cosmic star formation activity to present day 
observational evidence of their fossil remnants. 

A top-heavy IMF for the early episodes of star formation in the universe might have 
other interesting observational consequences, we speculate further on them below. The 
kinetic energy released during the thermonuclear explosions powered by pair instability are 
^10^ larger than those of ordinary Type II SNe. This might cause the interaction with the 
circumstellar medium to be as strong as predicted for hypernovae (Woosley & Weaver 1982). 
However, these explosions do not lead to the ejection of strongly relativistic matter (Fryer, 
Woosley & Heger 2001) and therefore cannot power a gamma-ray burst (GRB). 

In PopIII progenitors of VMBHs, the estimated angular momentum is sufficient to delay 
black hole formation and the system might develop triaxial deformations (Fryer, Woosley & 
Heger 2001). If the instabilities have enough time to grow, the core might break into smaller 
fragments that would then collapse and merge to form the central VMBH. If not, the star 
might still develop a bar-like configuration. Both these scenarios lead to a significant emission 
of gravitational waves (Schneider et al. 2000; Fryer, Holz & Hughes 2001). Furthermore, 
significant emission of gravitational radiation can occur as a result of the in-spiral and merger 
of VMBHs onto the SMBHs in the center of host systems (Madau & Rees, 2001) 

Once the VMBHs have formed, accretion continues through a disk at a rate which can 



-19- 



be as large as IOMqS"^ (Ftyer, Woosley & Heger 2001). Magnetic fields might drive an 
energetic jet which can produce a strong gamma-ray burst through the interaction with 
surrounding gas. The properties of these PopIII GRBs would be considerably different from 
their more recent {z < 5) counterparts: depending on the uncertain interaction of the jet 
with the surrounding matter, the bursts would be probably longer [10(1 + z) s] and the 
peak of emission, that in the rest frame is in 7-rays, would be shifted into X-rays. Indeed, 
BeppoSAX has revealed the existence of a new class of events, the so-called X-ray flashes or X- 
ray rich GRBs, which emit the bulk of their energy in X-rays (Piro, private communication). 
Furthermore, since PopIII GRBs explode at very high redshifts, it is hkely that their optical 
afterglow might be heavily absorbed by the intervening gas. These systems might be the 
natural candidates for the significant number (about 40 % of GRBs for which fast follow-up 
observations were carried out) of GRBs which do not show an optical counterpart, the so- 
called GHOST (GRB Hiding Optical Source Transient) or dark GRBs. Other explanations 
ascribe the failed optical detection to dust extinction within the host system but the ultimate 
nature of GHOSTS is still very debated (Lazzati, Govino & Ghisellini, 2001; Djorgovski et 
al. 2001) 

Finally, the energetic jets generated by GRBs engines produce, by photon-meson in- 
teraction, a burst of TeV neutrinos while propagating in the stellar envelope (Meszaros & 
Waxman 2001). We investigate this aspect in a companion paper (Schneider, Guetta & Fer- 
rara 2002) where we use the constraints set by the AMANDA-BIO experiment on the total 
integrated fiux of TeV neutrinos from PopIII GRBs. 
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Fig. 1. — Upper panel: The evolution of the temperature as a function of the hydrogen number 
density of proto-stellar clouds with the same initial gas temperature but varying metallicities Z = 
(0, 10~^, 10~*, 10~^, 1)^0 (Z increasing from top to bottom curves). The dashed lines correspond to the 
constant Jeans mass for spherical clumps; open circles indicate the points where fragmentation stops, filled 
circles mark the formation of hydrostatic cores. This Figure is reproduced from Figure 1 of Omukai (2000) 
for illustration after some modifications. Lower panel: The adiabatic index 7 as a function of the hydro- 
gen number density for the curves shown in the upper panel. Dotted (dashed) lines correspond to 7 = 1 
(7 = 4/3); open and filled circles as above. 
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Fig. 2. — Fractional abundance of H2 (solid lines) and ionization fractions (dashed lines) as a function of 
gas temperature. The three cases plotted are for a halo of mass M = IO^Mq at three different virialization 
redshifts z = 15,20,25, (virial temperatures of Tvir = [2.2,6.9,17] x 10^ K). As can be seen in the Figure, 
the evolution appears to be independent of the virial temperature of the halo. 
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Fig. 3. — Metal yields of the main elements synthesized in metal- free SN^^ according to Heger & Woosley 
(2001). The upper solid line corresponds to the total yield ejected and the dashed line indicates the kinetic 
energy of the explosion. 
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Fig. 4. — The top-heavy to normal IMF transition redshift, z/, as a function of SN-y^ progenitor mass 
fraction and the mass density contributed by VMBHs Qvmbh- Upper Panel: The computed critical density 
of VMBH remnants is compared to the observed value for SMBHs (upper dashed horizontal line). Lower 
Panel: The computed critical density of VMBH remnants is compared to the contribution to fl from the 
X-ray bright, off-center ROSAT sources (the lower dashed line) and to the abundance predicted assuming 
that the baryonic dark matter in galaxy halos is entirely contributed by VMBHs (upper dashed line). The 
observations on ^Ivmbh constrain the value for f^-^. For a given f^^, the corresponding value for the 
transition redshift can be inferred by the Zf curve. 




